Consistency in the driven butterfly 
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We study the consistency of the responses in the chaotic Lorenz model driven by chaotic and 
colored-noise signals. A measure is introduced for the consistency of the response in the system for 
repeated applications of the same drive signal, starting at different initial conditions of the response 
system. The convergence of responses is analyzed in dynamical phase space in terms of the local 
conditional Lyapunov exponent. The dependence of the local conditional Lyapunov exponent on 
drive signal, and the existence of multiple basins of consistency depending on strength of drive 
signal, are shown. 
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The dynamical response of complex structures rang- 
ing from biomolecules to the Millennium Bridge [l[ to 
complex drive waveforms is an important issue for sci- 
ence and engineering. Analysis of dynamical responses 
can provide a means for non-invasive testing and evalu- 
ation of the integrity of structures and changes due to 
damage or aging This approach may be extended 
to characterize the behavior of systems with complex in- 
ternal dynamics, such as electronic circuits, laser sys- 
tems, musical instruments and even brains. However, 
characterizing the responses of nonlinear dynamical sys- 
tems with complex dynamics presents many challenges. 
It may appear at first sight that the wide variety of possi- 
ble responses for a system with complex dynamics would 
make this task hopeless. On the other hand, a system as 
complex as the brain can have responses with various de- 
grees of reproducibility and consistency when subjected 
to repeated inputs of information and stimulation. It is 
our goal here to develop a new and practical approach 
to analyzing complex drive-response behavior based on 
measuring consistency |3[ of responses to repeated appli- 
cations of the same drive signal. 

Pecora and Carroll [H analyzed the stability of the re- 
sponse of a dynamical system to a drive signal through 
evaluation of the conditional Lyapunov exponent, condi- 
tioned on the drive signal. In their examples, the drive 
waveforms were generated by low dimensional chaotic 
system, and the drive and response systems displayed 
identical synchronization [f|. The more general case of 
stable synchronized response where drive and response 
waveforms are not identical has been referred to as gen- 
eralized synchronization Q . The condition of generalized 
synchronization corresponds to the condition for identical 
synchronization of two replica auxiliary systems driven 
by the same input signal. In the past decade, many 
studies have extended the drive signals to include noisy 
waveforms 0, [s| , which are very high dimensional. Re- 
cent research on coherence resonance has included exper- 
iments on colored- noise driven excitable systems Q . We 



also note that the concept of reliability has been used 
to describe the timing consistency of the evoked spike- 
sequences of rat cortical neurons in response to repeated 
presentations of the same noise waveform [loj . 

In this Letter, we consider the consistency viewpoint 
which emphasizes the reproducibility of responses to a 
particular input in driven nonlinear s yste ms. We use the 
specific example of the Lorenz system [Tl| driven by chaos 
and noise signals. We show the onset of consistency in 
terms of a consistency measure and the conditional Lya- 
punov exponent [J], which measures the convergence of 
trajectories in the space of responses conditioned on a 
particular drive signal. We also examine the local condi- 
tional Lyapunov exponent (LCLE) [l2| which allows us 
to see the contraction and expansion regions, and see the 
dependence of the distribution of the LCLE on changes 
of the drive signal. Finally, we consider the existence of 
multiple basins of consistent response for the same drive 
signal and the dependence of the basins on the drive sig- 
nal. 

Consider a dynamical system given by dX(t)/dt = 
F(X(t), S(t)), where X is the state variables of the re- 
sponse system, S is the drive signal, and F is a nonlinear 
function. We consider one of the observable state vari- 
ables x(t). The consistency C can be measured as the av- 
erage correlation between responses for a set of response 
system initial conditions, 
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where x(f;Xj(0)) and x(t;Xj(0)) are variables x(t) start- 
ing from different initial conditions of Xj(0) and Xj(0), 
Xi and Xj are the mean values of the two time series of 
x(i;Xj(0)) and x(t;Xj-(0)), and <Ji and <7j are the stan- 
dard deviations of the two time series. The angle brackets 
denote time averaging. 
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We use the Lorenz model ll| with an additive drive 
for the y- variable, dy/dt = —xz + rx — y + Ds(t), where 
the drive signal s(t) is multiplied by the coefficient of 
drive strength D. Two examples of s(t) are (i) y variable 
of the Lorenz chaos generated from the basic model and 
(ii) an exponentially correlated colored noise waveform 
generated using an Ornstein-Uhlenbeck process [l3| with 
a correlation time r c = 0.2. The temporal waveforms of a 
drive signal and corresponding response of the x variable 
starting from ten different initial conditions are shown in 
Figs. 1(a) and 1(b) for the chaos and colored-noise drive 
signals. It can be seen that the waveforms of the re- 
sponses converge after a short transient. We calculated 
the consistency measure C of Eq. (TTJ) for an ensemble 
of initial states randomly selected on the unperturbed 
Lorenz butterfly. Figures 1(c) and 1(d) show the asymp- 
totic value of C (— » 1) as a function of the drive strength 
D for the Lorenz chaos and colored-noise drive signals. 

To calculate the conditional Lyapunov exponent we 
begin with the linearized equations of d£(t)/dt — 
D x F(£(i),X(t)), where £ are the linearized variables of 
X and D X F is the Jacobian matrix with respect to X 
(not including S). The conditional Lyapunov exponent 
A c for consistency is defined as, 
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where ||£|| is the norm of the linearized variables. We only 
consider the maximal value of the conditional Lyapunov 
exponent for simplicity. We calculated A c and plotted 
it as a function of the drive strength D in Figs. 1(c) 
and 1(d). It can be seen that the onset of consistency 
[C — 1) corresponds to a negative value of the max- 
imum conditional Lyapunov exponent, which indicates 
that close trajectories tend to converge on the average 
when the system is driven strongly enough. 

For comparison, we also calculated the usual Lyapunov 
exponent A by using the linearized equations of the whole 
drive-response system for the chaos-driven Lorenz model; 
dr)(t)/dt = D x , s F(»7(t),X(t),S(i)), where rj are the lin- 
earized variables of both X and S, and D x S F is the Ja- 
cobian matrix with respect to both X and S. A cannot 
be calculated for noise-driven systems because of infinite 
dimensionality of the noise drive signal S(t). A is plotted 
as a function of the drive strength D as shown in Fig. 
1(c). It is worth noting that the measure of consistency 
A c is different from the measure of chaoticity A in the 
presence of drive signals. At large drive strengths consis- 
tent behavior can be observed (A c < 0) even for a chaotic 
motion (A > 0). 

In order to see the convergence behavior locally in time 
and space, we introduce the local conditional Lyapunov 
exponent (LCLE), 
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where At is a finite time corresponding to the time step 
of the numerical calculation. £(f) is calculated along the 



trajectory of the response system, which is influenced by 
the time-dependent drive signal. ||£(i)|| is normalized 
at each step of numerical integration so that the vector 
£(t) stays a unit vector. It is important to maintain the 
direction of the vector £(t) to incorporate the influence 
of the time-dependent trajectory in the phase space. 

We plotted A; c on the trajectory of the Lorenz 'butter- 
fly' attractor on the x-z plane of the phase space. Figure 
2(a) shows A; c before a drive signal is added. Figures 2(b) 
and 2(c) show A; c for a response driven by a Lorenz chaos 
and a colored noise drive signal, respectively. Strong 
expansion (A; c > 0) and contraction (A; c < 0) regions, 
are indicated by red and blue regions. The red regions 
show the locations on the attractor where the system 
is more easily susceptible to breakdown of consistency 
due to perturbative noise. Comparing the Figures 2(a), 
(b) and (c) it can be seen how the distribution of con- 
verging and diverging dynamics can depend on the drive 
signal. Zhou and Kurths discussed contraction regions in 
the driven butterfly in the context of synchronization 
using a definition of contraction region based on he Ja- 
cobian matrix of the response system, rather than the 
LCLE. This approach does not incorporate the depen- 
dence of the convergence on the time dependent nature 
of the drive signal. 

We show the change in shape of the distribution of Aj c 
values for increasing drive strength D in Fig. 3(a). The 
distributions of A/ c at constant drive strengths D = 0, 
60 and 120 are shown in Figs. 3(b)-3(d). These figures 
correspond to the average values of A c = 0.905, -0.110 
and -1.507, respectively. The shape of the distribution 
shifts towards negative values of A; c with increase of drive 
strength. However, the position of the main peak, cor- 
responding to weak contractions, does not change, and 
the onset of consistency does not require a drastic shift 
of the whole distribution to the negative side. 

As a final point of discussion, we note that the LCLE 
measures the convergence of close trajectories, and does 
not guarantee the convergence of trajectories from arbi- 
trary positions in the phase space. A response system 
could have multiple basins - that could be more than one 
response attractor. An example is the case of r = 13, 
where the Lorenz system has two fixed points without 
a drive signal. When we drive the x variable with a 
colored-noise signal of small drive strength, D = 10, (and 
r c = 2.0) we can obtain two types of response trajecto- 
ries for the same drive signal depending on the initial 
conditions. Figure 4(a) shows two separate response tra- 
jectories (black curves) located near the original fixed 
points on the x-y plane of the phase space. The color 
code shows the probability of converging into the con- 
sistent trajectory in one of the two basins for a set of 
colored-noise drive signals (ten different temporal wave- 
forms). The basins of convergence for each of the two 
possible responses are shown in red and blue. Each of 
these responses can be seen to be locally consistent in 
the sense that trajectories starting in the same basin con- 
verge after transients. The boundary of the two basins 
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depends on the drive signal. This example shows that 
the ensemble of initial states is significant for the evalua- 
tion of consistency measures as defined in Eq. (1). When 
the drive strength is increased to D = 30 the amplitudes 
of the two attracting response trajectories become larger 
and merge together. For large drive strengths there ex- 
ists only one attracting response trajectory and a single 
basin of consistency for all initial conditions, as shown in 
Fig. 4(b). 

In conclusion, we have introduced the viewpoint of con- 
sistency in driven nonlinear dynamical systems and used 
the Lorenz model driven by chaos and colored noise sig- 
nals to illustrate this viewpoint. Consistency is evalu- 
ated by the cross-correlation function between responses 
starting from different initial conditions when the system 
is repeatedly driven by the same waveform. Consistent 
responses are observed after transients even though the 
system is in a chaotic state. A rather different insight 
is provided when we compute and plot the LCLE in the 
phase space, which allows us to visualize contraction and 
expansion regions on the driven chaotic attractors. We 



observe the dependence of the distribution of the LCLE 
on changes of the drive signal to understand the onset of 
consistency. We also found the multiple basins of consis- 
tency, where the Lorenz system exhibits more than one 
type of response with respect to the same drive signal, 
and yet is still consistent in the local sense that the sys- 
tem converges to the same response for large sets of initial 
states. We expect that the viewpoint of consistency in 
driven nonlinear dynamical systems will be useful for an- 
alyzing and understanding drive-response characteristics 
of complex systems in science and engineering. 
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FIG. 1: (a),(b) Temporal waveforms of the drive signal and 
the x variable of the driven Lorenz model starting from ten 
different initial conditions, (a) Chaos drive and (b) colored 
noise drive. (c),(d) The consistency C (black curve), the con- 
ditional Lyapunov exponent A c (blue curve), and the Lya- 
punov exponent A (red line) as a function of drive strength 
D. (c) Chaos drive and (d) colored noise drive, a — 10, 
r = 28, and b — 8/3 are often used for the Lorenz model to 
observe chaotic behaviors without a drive signal, y variable 
of the Lorenz model is additively driven by a Lorenz chaos 
and a colored noise signal with a correlation time r c = 0.2. 



FIG. 2: The driven butterfly, (a) No drive, (b) chaos drive 
(D = 15), and (c) colored noise drive (D = 70). The color 
indicates the local conditional Lyapunov exponent plotted on 
the trajectory in the x-z plane of the phase space. Expan- 
sion (red) and contraction (blue) regions are observed on the 
trajectory. 
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FIG. 3: (a) Three dimensional picture of the distribution 
of local conditional Lyapunov exponent Xi c as a function of 
colored-noise drive strength D. (b)-(d) Distribution of A; c at 
constant drive strengths, (b) D = 0, (c) D = 60, and (d) 
D = 120. The conditional Lyapunov exponent A c is obtained 
from the average of the distribution of A; c : (b) A c = 0.905, 
(c) A c = -0.110, and (d) A c = -1.507. 



FIG. 4: (a) Multiple basins of consistency at D = 10. Black 
lines show attracting trajectories which correspond to locally 
consistent responses. Color indicates the probability of con- 
verging into a trajectory in the basin 1, plotted on the x — y 
plane of the three-dimensional phase space at z = 12. Red 
corresponds to the basin 1 and blue corresponds to the basin 
2. (b) Single basin of consistency at D = 30. Black line in- 
dicates the consistent trajectory. Red indicates the basin 1. 
The Lorenz model has two fixed points without a drive signal 
(r = 13). a;- variable of the Lorenz model is additively driven 
by a colored-noise signal with a correlation time r c = 2.0. 
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